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The spontaneous breaking of chiral symmetry (SBCS) is one of the most important phenomena 
of hadron physics. It defines the properties of all the light mesons and baryons. The Chiral Per- 
Q_i turbation Theory (ChPT) encodes QCD hadronic correlators at low-energy region in the terms 

O | the low-energy constants (LEC) - the expansion parameters on light quark current masses m and 

external momenta p. The LEC's can be extracted from the phenomenology or from QCD lat- 
tice calculations. On the other hand, QCD instanton vacuum/instanton liquid model provides a 
very natural nonperturbative explanation of the SBCS. It provides a consistent framework for de- 
scription of the pions and thus may be used for evaluation of the LEC. Our aim is to calculate the 
vacuum properties and the LEC's within instanton vacuum model and confront with phenomenol- 
ogy and lattice results. 

o 

(N 



43 



> 

X 



XXI International Baldin Seminar on High Energy Physics Problems, 
September 10-15, 2012 
JINR, Dubna, Russia 



* Speaker. 

© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-Share Alike Licence. http://pos.sissa.it/ 



IR QCD and SBCS 



Mirzayusuf Musakhanov 



1. Introduction 

1.1 Nonperturbative QCD in Infrared Region 

QCD in the chiral limit is a good approximation to the real world and left and right-handed 
quarks are decoupled. But the hadrons has no chiral doublets which means that the QCD vacuum 
breaks the chiral symmetry. One of the signals of the Spontaneous Breaking of the Chiral Symme- 
try (SBCS) is the presence of the nonzero chiral quark condensate, (qq) / 0. The understanding of 
SBCS is provided by Casher-Banks formula 1 : 

(??) = ~v(A = 0), (1.1) 

The chiral condensate is thus proportional to the averaged spectral density of the QCD Dirac op- 
erator (id + gA) at zero eigenvalues v(A = 0) [1, 2]. The key moment is that the Dirac operator 
in the background of topologically nontrivial field may has an exact zero modes with A = and 
in an accordance with a general Atiah-Singer index theorem the number of these modes equal to 
Pontryagin index or the topological charge of the background field. Then, the SBCS is due to of a 
derealization of the would-be zero modes, induced by the background of a vacuum topologically 
nontrivial fields, resulting from quarks hopping between them. On the other hand the QCD sum 
rules phenomenology requests that the QCD vacuum has the gluon condensate [3] : 

3^2<G£ V G£ V )~ (200 MeVf. (1.2) 

The simplest way to explain both phenomena - the chiral quark condensate and the gluon conden- 
sate is the instanton vacuum/instanton liquid model (see reviews [1, 2]). In the model the gluon 
condensate is directly related to the instanton density as 3^1 (G^ V G^ V ) = N/V = R , where from 
Eq.(1.2) the average inter-instanton distance in Euclidian space R ~ (200MeV)~ l = 1/m. With 
this value for R and phenomenological value for the quark condensate (qq) ~ — (250 MeV) 3 the 
average size of the instanton is estimated as p ~ 0.33 fm. 

The main problem of the instanton vacuum model is the lack of the confinement. A gen- 
eralizations of the model providing the confinement was proposed recently with the price of the 
introducing of the other topological objects. Among them are calorons with non-trivial holonomy, 
BPS monopoles or dyons, vortices, etc.. 

1.2 Lattice QCD in Infrared Region 

In the recent decades QCD vacuum have been intensively studied by direct numerical simu- 
lations on the lattice. The presence of non-trivial topological objects was demonstrated by using 
various configuration-smoothing methods. The typical view of the distributions of the action and 
topological charge density before and after few steps of the "cooling" was the following: before 
"cooling" they are heavily dominated by perturbative zero-point fluctuations of the gluon fields. 
The "cooling" suppresses these fluctuations which leads to a smooth background coinciding with 
an ensemble of instantons and antiinstantons with p ~ 0.36 fm and R ~ 0.89/m [4]. 

'Assumed thermodynamic limit means that the volume V goes to infinity faster than quark mass m goes to zero. 
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But recently the QCD instanton vacuum model was challenged by the detailed lattice studies of 
topological charge density distribution by means of overlap Dirac operator (possess an exact chiral 
symmetry on the lattice) which shows a three-dimensional, laminar and highly singular structure 
which seems contradict the instanton and similar pictures [5,6]. Probably the structure of the QCD 
vacuum much more rich and we have to add another topological objects mentioned above. 

The question: which component of QCD vacuum is most important for the SBCS? The answer 
is given by the comparison of the momentum dependence of the dynamical quark mass M{q) in 
the chiral limit given at instanton vacuum model and lattice QCD, Fig.l. Model ^-dependence of 
M(q) is mostly due to the quark zero-mode function in the (anti)instanton field. We see very good 
coincidence without any fitting of the lattice results with the model. 

OA 




Figure 1: Momentum dependence of dynamical quark mass M(q) in the chiral limit. Points: lattice result 
[7]. Red line: instanton vacuum model [8], no fitting. 

The details of the instanton size distribution extracted from lattice simulations of the instanton 
liquid are presented at the Fig. 2. 
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Figure 2: » mode /(p) « exp 
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with p = 0.3 fm, a = 0. 1 3 fm, R ~ 1 .07 fm and the VMP distribution 



from the corresponding lattice configurations [9]. 
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1.3 QCD vacuum and pion physics 

The properties of the QCD vacuum and its lowest excitations-pions are encoded in Low- 
Energy Constants (LEC) of Chiral Perturbation Theory (ChPT). It is natural to confront LEC's 
derived from the model with the phenomenology and lattice calculations of them. ChPT gives the 
correlators in terms of LEC's as coefficients in the expansion over the current quark mass m and 
external momentum of the order of M K (taking into account that in the lowest order m\ ~ ra). This 
expansion have to take into account a chiral log terms ~ log ra 2 due to the contribution of the pion 
loops. Nf = 2 lowest q 2 order effective lagrangian L2 have two LEC's F and B. They correspond 
to the observables: F - the pion decay constant in the chiral limit and ra 2 = 2Bm - the pion mass 
square in the lowest order on ra [10]. At q order effective lagrangian L4 has 10 independent bare 
LEC's lj,hj. They are renormalized by pion loops to Tj. Physical observables should be expressed 
in terms of T\ [10]. 

Our aim is to calculate the vacuum properties and the LEC's within instanton vacuum model 
and compare with phenomenology and lattice results. 



2. Instanton vacuum model 

Instanton A 1 (antiinstanton A 1 ) is a solution of Yang-Mills equations in Euclidian space with 
topological charge = 1(— 1) and correspond to the tunneling process of the gluon fields. Instanton 
collective coordinates are 4 (centre) + 1 (size) + (4N C — 5) (orientations) = 4N C . It is assumed 
that vacuum background is given by A = £/A 7 + L/A 7 . The main parameters of the instanton 
vacuum are average interinstanton distance R and average size p. The estimates of them are the 
following: lattice estimates: R ps 0.89/ra, p ps 0.36 fm [4], phenomenological estimates: R ps 1/m, 
p ps 0.33 fm [1], our estimate (corresponding ChPT F n , m= Q = 88MeV, (qq) m=0 = -{255MeV) 3 ): 
R ~ 0.76/m, p « 0.32 fm, Thus within 10 — 15% uncertainty different approaches give similar 
estimates. We see that the packing parameter 7T 2 (|) 4 ~ 0.1 is small and independent averaging 
over instanton positions and orientations is justified. 

2.1 Light quarks in the instanton vacuum 

We have to calculate the correlators beyond the chiral limit to extract LEC's. Our starting 
point is the interpolation formula [11, 12, 13, 14, 15] for the quark propagator in the field of single 
instanton: 

St = So-Sop J^?,°i x PSo, So = t-Jt- , (id +gAi)<i> 0i = (2.1) 
(9oi\pS p\9oi) p + im 

where <J>o, is the quark zero-mode function. The advantage of the formula shown by the projection 
of the propagator to the zero-mode 

Si\*0i) = —\^oi), (®0i\Si = . (2.2) 

im im 

as it must be. Also, the calculations of the single-instanton quark effective action reproduce exactly 
most important for our problems the (rap ) 2 In mp -term from [16]. 
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On the other hand the interpolation formula provide the summing up the re-scattering series 
for the full propagator (in the presence also of the flavor external fields V = s + py$ + v + ay*,): 



S-S = -So£p\<l>Qi)($Oi (^rr) 



<foj ) (fojlPSo, (2.3) 



\(j> ) = -^Lp\®o), §0= } , Li(x,Zi) =Pexp( i f <fy M (v M (y) + a^(y)y 5 )) . 
p p + V + im V Jzi J 

Gauge connections L, provide flavor gauge-covariance of the propagator S. 

Let's represent the quark determinant as Det = Det/^ x Det/ ow [8]. The high-frequencies part 
of the quark determinant Det/,^ can be calculated in single-instanton approximation, while the 
low-frequencies part Det/ mv . get the contribution from the whole set of the instantons according the 
formula: 

lnDet w = Tr / dm' S(m') = lndet{^,i\pS f 8 p\^ QJ ) , (2.4) 

J m 

The partition function for the light quarks Zn [V] is given by the averaging of Det/ 0lv over instantons 
collective coordinates. It was done by means of the fermionization with constituent quarks yf , ¥ 
and leads to the t'Hooft-like nonlocal quark interaction term J#, with 2Nf legs [11, 12]. In the 
following we consider the Nf = 2 case, neglecting the influence of strange and heavier quarks for 
the pion physics observables. 

Further exponentiation in Zn[V] introduce the integration over dynamical couplings A± and 
we have the partition function in the form: 

Z N [V] = J dX+dX-D^DYe^, (2.5) 

S=xi/^id+V + im)y+Y,\N ± la^--N ± + X±yA , 



2N r -\ „o a 



2 



(N}-\)2Nc 8A^ c -4' 
Jf g = (2np) 2 ¥ }LF ] ^FLxtr g , J± vfg = {Inpf^LFa^ ^FL %J L = y^y 4 . 



The interaction term non-locality form-factor F in the momentum space is completely defined by 
Fourier-transform of the quark zero-mode function as: 

F ( k)=- d j m m)-h(>)Km,^ I ^ m ,np<3-, (j^,i*iP>3. (2.« 

So, the range of the non-locality is given by the instanton size p, as was expected. 

The main purpose of our work is the calculations of various correlators with account of the 
ff(l/N c , m, m/N c , m/N c lnm)-corrections, which means double expansion over l/N c and m. We 
estimated suchNLO corrections [13, 14, 15, 17]: 

1. The width of the instanton size distribution is G{\ //V c ).We found that the finite width correc- 
tions are negligible, as were expected. For example, the account of the finite width parameter 
a (see n mo d e i{p) from the Fig. 2) lead to the corrections « 2.6% for (qq) and « 5% for F%. 
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2. The back-reaction of the light quark determinants to the instanton vacuum properties does 
not sizably change the distribution over N + + N- but radically change the distribution over 
N + —N-. Any my = leads to c) -function type of the latter distribution. In the following we 
take N+ = AL_. 

3. There are the quark-quark tensor interaction terms which are 1/A^-suppressed and thus are 
absent in the LO. It might give the contribution to the vacuum magnetic susceptibility %. Our 
estimate shows that it leads to the correction < 1% to the %, which is certainly negligible. 

4. We found the importance of the meson loops contribution. And most important among them 
are pion loops, certainly. They lead to the chiral logs and numerically large corrections. 

Though the evaluation of the meson loop corrections in the instanton vacuum model is similar to 
the earlier meson loop evaluations in the NJL model [18, 19, 20, 21] important differences should 
be mentioned: 

1. Due to nonlocal formfactors there is no need to introduce independent fermion and boson 
cutoffs A/, Aft. The natural cutoff scale for all the loops (including meson loops) is the 
inverse instanton size p ~ 1 . 

2. In the instanton vacuum model the quark coupling constant X is defined through the saddle- 
point equation whereas it is a fixed external parameter in NJL. 

3. Vacuum properties 

First it were calculated the dependencies of the dynamical quark mass M and the quark condensate 
(qq) on current quark mass m [13], see Fig. 3. 



-<99>, GeV 3 

M(m), GeV 




0.025 0.05 0.075 0.1 0.025 0.05 0.075 0.1 

Figure 3: Left panel: m-dependence of the dynamical quark mass M on the scale p _1 sa Q.6GeV. The 
solid curve - the exact numerical solution Eq. (3.1). The dashed curve - the solution obtained by the 1 /N c 
iterations with the same accuracy. Lattice data points are from [22], where the scale is 1.64 GeV. 
Right panel: ra-dependence of the quark condensate — (qq) Eq. (3.2). The long-dashed curve is the LO 
result, the short-dashed curve is the NLO contribution, the solid curve is the total one. The dot-dashed line 
is the leading-order in 1 /A^ c -expansion. 
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We found M(m) (in GeV) and (qq)(m) (in GeV ) as a functions of m (in GeV) as 

M(m) = 0.36-2. 36m-— (0.808 + 4.1971nm) + <^ (m 2 , — ] ,(3.1) 

-(qq) (m) = (0.00497 - 0.034m) N c + (0.00168 - 0.049m - 0.058m lnm) + ff (m 2 ,^j (3.2) 

We see a very essential contributions of the chiral logs to the m dependencies of M(m) and (qq) (m). 
They are coming from the pion loops, certainly. 

3.1 Vacuum magnetic susceptibility 

External electromagnetic field F^ v generate the quark currents in the QCD vacuum, described 
by the magnetic susceptibility %f( m f) (with normalization factor (i\jf^\jf)o = —(qq) (m = 0)): 

(0|y^ct M vV/|0>f = e f Zf{m f ) (z>V}o^>v (3-3) 

QCD sum rules [23, 24] gave Xf( m f) (W^W)o ~ 40 — 70MeV, while recent lattice measurements 
gave for this quantity « 46(3) MeV (quenched SU(2) QCD, chiral limit) [25], « 52MeV (quenched 
SU(3) QCD, chiral limit) [26], for (u,d)-quarks 40(1 .4) MeV and for s-quark » 53 (7.2) MeV 
(fully dynamical QCD) [27]. On the other hand, the magnetic susceptibility is measurable in jet 
production at q± 3> Aqcd at the process y+N — > (qq) +N with polarized photon [28]. 
Our result for %(m) (iyfixfio [15] with the accuracy 6 [n? , ^ J is 

X(m) (iY f ¥)o = N c (o.015 + 5 • 10~ 4 m + ^ lnm ) ~ 0007 ~ 0.415ra - 0.198m lnm. (3.4) 



X<«>, GeV 



i<vi) NW 



I 1 1 1 1 m, GeV 

0.025 0.05 0.075 0.1 



Figure 4: Magnetic susceptibility %{m) (i\j/< \jf)o as a function of current quark mass m. 
Last term in (3.4) is in the correspondence with the newly established chiral log theorem: 

,2 



XH = X(0) (l - ^272 ln "^) (3-5) 
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We see that the strange quark contribution to the magnetic susceptibility is much less than lightest 
quarks one, Fig. 4. This ra-dependence contradict to the ra-dependence obtained at [27] and the 
problem needs further research. 



4. Pion properties and low-energy constants I3, 14 

We calculated the two-point axial-isovector currents correlator [15]: 



q^qv 



+ 



(4.1) 



Here M n has a meaning of pion mass and F n - pion decay constant. 
We found 



F*=Nc \ I 2.85 



0.869 
1.63 



3.51 



0.815\ 44.25 



m 



^ = -((3-49 + — +m 15.5 + — + 



N c J N c 
18.25 13.5577 



mlnm + % 2 ) • 10~ J [GeV 1 ] 



lnm + G{m 2 ) [GeV 2 ] (4.2) 



represented by Fig. 5. 



F GeV 




0.075 0.1 




0.025 0.05 



0.075 0.1 



Figure 5: Left panel: m-dependence of the pion decay constant F n . The long-dashed curve is the LO con- 
tribution, the short-dashed curve is the NLO contribution, the solid curve is the total LO+NLO contribution. 
The dot-dashed line represents the leading-order in 1 /A^ -expansion result, evaluated with the mass Mo. 
Right panel: m-dependence of the pion mass M K . The long-dashed curve is the LO contribution, the short- 
dashed curve is the NLO contribution, the solid curve is the total LO+NLO contribution. The dot-dashed 
line represents the leading-order in 1 /Af c -expansion result, evaluated with the mass Mq. 



According to [10], the low-energy constants I3, I4 of the chiral lagrangian may be extracted 
from the <??(ra) -corrections to physical quantities, e.g. 



M% = ml I 1 



in- 



1 + 



//;; 



SK 2 F 2 



(4.3) 



Here m\ = 2mB. Lowest order LEC's B = 2.0 19 GeV and F = SSMeV were taken as an input to 
fix main instanton vacuum parameters as p = 0.350/m and R = 0.856/m. 
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Then, we found (neglecting by the terms & \ w)^ 

l 3 = -1 . U25N C + 0.0738 - 0.999 lnm, f 4 = -0.0793N C + 0.01876 - lnm, 



(4.4) 



which gives at m = 0.0055 GeV, h = 1.84, l 4 = 4.98 and corresponding M n = 0.142 GeV, F n = 
0.0937 GeV. Our values of (I3J4) should be compared with the phenomenological estimates [10, 
29] as well as lattice predictions [30, 31, 32, 33, 34, 35, 36] given in Fig. 6. 



JLQCD/TWQCD 2008, N f = 2 

RBC/UKQCD 2008 

PACS-CS (preliminary) 

JLQCD 2007, N, = 2 

ETM 2007, N, = 2 

Del Debbio et al. 2006, N, = 2 

MILC 2007 

Gasser & Leutwyler 1 984 
Instanton vacuum model 



JLQCD/TWQCD 2008, N, = 2 
RBC/UKQCD 2008 
PACS-CS (preliminary) 
JLQCD 2007, N, = 2 
ETM 2007, N, = 2 
MILC 2007 

Colangelo, Gasser, Leutwyler 
Gasser & Leutwyler 1 984 
Instanton vacuum model 



Figure 6: The LEC's It, (left panel) and I4 (right panel) - recent lattice results from different collabora- 
tions [30, 31, 32, 33, 34, 35, 36], phenomenological estimates [10, 29] and our result. 



4.1 Isospin breaking by m u — / and LEC's hi, and h 

The correlator of the iso-vector pseudoscalar and iso-singlet pseudoscalar currents and split- 
ting between the (uu) and (dd) quark condensates in QCD are proportional to m u —nid and de- 
scribed by LEC's I13 and lj as 

i f dx(P 3 (x)P { \0))e iqx = %B *^~™d)h +0 fa2) ; ( Uu ) _ ^ = 4B 2 (m u -m d )h 3 . (4.5) 
LEC h define the pion masses shift due to m u — 7^ 0: 

9 z?2 

(M 2 KQ -Ml ± ) Sm = -(m^mtf—l-j « -{m u - m d f\2 ■ W%. (4.6) 

Phenomenological estimate (GL, AP84) l-j ~ 5 ■ 10~ 3 . Our estimates [17] for 1-j and are 

/ 7 «(6.6±2.4)-10" 4 , /i 3 ~ 5.48 - 10 3 . (4.7) 

We found strong dependence of the lj on the instanton vacuum parameters p and R. Confronting 
with the possible lattice calculation of this one will fix these parameters well. 



5. Summary 

SBCS is generated by the 'hopping' of quarks between topologically nontrivial gluon lumps 
leading to changing of quark chirality. The strong evidence in favor of the instanton vacuum model 
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given by the momentum dependence of the dynamical quark mass M(p) found in the lattice calcu- 
lations (see Fig.l). 

The proposed interpolation formula (2. 1) for the quark propagator in the single instanton field 
and in the presence of the flavor external fields lead to the low-frequencies quark determinant (2.4). 

Next task was a derivation of the light quarks partition function Z[V], which is a functional of 
the external flavor fields V. For this one we averaged the quark determinant over instantons, using 
diluteness of the instanton media. On this way it was introduced constituent quarks as a tool of this 
averaging. 

The partition function (2.5) has t'Hooft-like interaction terms with 2Nf quark legs, where the 
coupling X is not a fixed one, but must be calculated from saddle-point condition. We considered 
Nf = 2-case. 

Since we are going beyond of the chiral limit we have to calculate the correlators with account 
of 0(l/N c ,m,m/N c ,m/N c \nm)-conections. Within the model there are the contributions of meson 
loops, finite width of instanton size distribution and quark-quark tensor interactions term to these 
corrections. We found that the most important are meson loops, while other contributions are 
negligible. 

ChPT describes SBCS in low-energy region in the terms of LEC's. We fixed the main param- 
eters of the model R and p in terms of lowest p 2 order LEC's B and F. With these parameters 
we also considered the vacuum magnetic susceptibility %f an d its dependence on the current quark 
mass nif, see Fig. 4. It might be important for the phenomenology of the jet production at the 
process y+ N — > (qq) +N with polarized photon. 

At further step we calculated next p 4 order LEC's J3, I4, I7 and h^. Most important among 
them are Z3 and I4, since there are available phenomenological and lattice estimates of them. We 
confronted our model calculation with available data and found reasonable correspondence, see 
Fig. 5. 

This means that the instanton vacuum is applicable for understanding of the low-energy hadron 
physics, at least on the qualitative level. 

I would like to thank Prof. Ernst-Michael Ilgenfritz for useful discussions during ISHEPP2012. 
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